LOSCAR: Long-term Ocean-atmosphere-Sediment CArbon cycle Reservoir Model v2.0.4

The LOSCAR model is designed to efficiently compute the partitioning of carbon between ocean, atmosphere, and sediments on time scales ranging from centuries to millions of years. While a variety of computationally inexpensive carbon cycle models are already available, many are missing a critical sediment component, which is indispensable for long-term integrations. One of LOSCAR's strengths is the coupling of ocean-atmosphere routines to a computationally efficient sediment module. This allows, for instance, adequate computation of CaCO<sub>3</sub> dissolution, calcite compensation, and long-term carbon cycle fluxes, including weathering of carbonate and silicate rocks. The ocean component includes various biogeochemical tracers such as total carbon, alkalinity, phosphate, oxygen, and stable carbon isotopes. LOSCAR's configuration of ocean geometry is flexible and allows for easy switching between modern and paleo-versions. We have previously published applications of the model tackling future projections of ocean chemistry and weathering, <i>p</i>CO<sub>2</sub> sensitivity to carbon cycle perturbations throughout the Cenozoic, and carbon/calcium cycling during the Paleocene-Eocene Thermal Maximum. The focus of the present contribution is the detailed description of the model including numerical architecture, processes and parameterizations, tuning, and examples of input and output. Typical CPU integration times of LOSCAR are of order seconds for several thousand model years on current standard desktop machines. The LOSCAR source code in C can be obtained from the author by sending a request to [email protected]

CO2 perturbation experiments: similarities and differences between dissolved inorganic carbon and total alkalinity manipulations

Increasing atmospheric carbon dioxide (CO2) through human activities and invasion of anthropogenic CO2 into the surface ocean alters the seawater carbonate chemistry, increasing CO2 and bicarbonate (HCO3−) at the expense of carbonate ion (CO32−) concentrations. This redistribution in the dissolved inorganic carbon (DIC) pool decreases pH and carbonate saturation state (Ω). Several components of the carbonate system are considered potential key variables influencing for instance calcium carbonate precipitation in marine calcifiers such as coccolithophores, foraminifera, corals, mollusks and echinoderms. Unravelling the sensitivities of marine organisms and ecosystems to CO2 induced ocean acidification (OA) requires well-controlled experimental setups and accurate carbonate system manipulations. Here we describe and analyse the chemical changes involved in the two basic approaches for carbonate chemistry manipulation, i.e. changing DIC at constant total alkalinity (TA) and changing TA at constant DIC. Furthermore, we briefly introduce several methods to experimentally manipulate DIC and TA. Finally, we examine responses obtained with both approaches using published results for the coccolithophore Emiliania huxleyi. We conclude that under most experimental conditions in the context of ocean acidification DIC and TA manipulations yield similar changes in all parameters of the carbonate system, which implies direct comparability of data obtained with the two basic approaches for CO2 perturbation

Early detection of ocean acidification effects on marine calcification

Ocean acidification is likely to impact calcification rates in many pelagic organisms, which may in turn cause significant changes in marine ecosystem structure. We examine effects of changes in marine CaCO3 production on total alkalinity (TA) in the ocean using the global biogeochemical ocean model HAMOCC. We test a variety of future calcification scenarios because experimental studies with different organisms have revealed a wide range of calcification sensitivities to CaCO3 saturation state. The model integrations start at a preindustrial steady state in the year 1800 and run until the year 2300 forced with anthropogenic CO2 emissions. Calculated trends in TA are evaluated taking into account the natural variability in ocean carbonate chemistry, as derived from repeat hydrographic transects. We conclude that the data currently available does not allow discerning significant trends in TA due to changes in pelagic calcification caused by ocean acidification. Given different calcification scenarios, our model calculations indicate that the TA increase over time will start being detectable by the year 2040, increasing by 5-30 mu mol kg(-1) compared to the present-day values. In a scenario of extreme reductions in calcification, large TA changes relative to preindustrial conditions would have occurred at present, which we consider very unlikely. However, the time interval of reliable TA observations is too short to disregard this scenario. The largest increase in surface ocean TA is predicted for the tropical and subtropical regions. In order to monitor and quantify possible early signs of acidification effects, we suggest to specifically target those regions during future ocean chemistry surveys

Marine biogeochemical responses to the North Atlantic Oscillation in a coupled climate model

In this study a coupled ocean-atmosphere model containing interactive marine biogeochemistry is used to analyze interannual, lagged, and decadal marine biogeochemical responses to the North Atlantic Oscillation (NAO), the dominant mode of North Atlantic atmospheric variability. The coupled model adequately reproduces present-day climatologies and NAO atmospheric variability. It is shown that marine biogeochemical responses to the NAO are governed by different mechanisms according to the time scale considered. On interannual time scales, local changes in vertical mixing, caused by modifications in air-sea heat, freshwater, and momentum fluxes, are most relevant in influencing phytoplankton growth through light and nutrient limitation mechanisms. At subpolar latitudes, deeper mixing occurring during positive NAO winters causes a slight decrease in late winter chlorophyll concentration due to light limitation and a 10%–20% increase in spring chlorophyll concentration due to higher nutrient availability. The lagged response of physical and biogeochemical properties to a high NAO winter shows some memory in the following 2 years. In particular, subsurface nutrient anomalies generated by local changes in mixing near the American coast are advected along the North Atlantic Current, where they are suggested to affect downstream chlorophyll concentration with 1 year lag. On decadal time scales, local and remote mechanisms act contemporaneously in shaping the decadal biogeochemical response to the NAO. The slow circulation adjustment, in response to NAO wind stress curl anomalies, causes a basin redistribution of heat, freshwater, and biogeochemical properties which, in turn, modifies the spatial structure of the subpolar chlorophyll bloom

Complementary constraints from carbon (13C) and nitrogen (15N) isotopes on the glacial ocean's soft-tissue biological pump

A three-dimensional, process-based model of the ocean's carbon and nitrogen cycles, including 13C and 15N isotopes, is used to explore effects of idealized changes in the soft-tissue biological pump. Results are presented from one preindustrial control run (piCtrl) and six simulations of the Last Glacial Maximum (LGM) with increasing values of the spatially constant maximum phytoplankton growth rate μmax, which accelerates biological nutrient utilization mimicking iron fertilization. The default LGM simulation, without increasing μmax and with a shallower and weaker Atlantic Meridional Overturning Circulation and increased sea ice cover, leads to 280 Pg more respired organic carbon (Corg) storage in the deep ocean with respect to piCtrl. Dissolved oxygen concentrations in the colder glacial thermocline increase, which reduces water column denitrification and, with delay, nitrogen fixation, thus increasing the ocean's fixed nitrogen inventory and decreasing δ15NNO3 almost everywhere. This simulation already fits sediment reconstructions of carbon and nitrogen isotopes relatively well, but it overestimates deep ocean δ13CDIC and underestimates δ15NNO3 at high latitudes. Increasing μmax enhances Corg and lowers deep ocean δ13CDIC, improving the agreement with sediment data. In the model's Antarctic and North Pacific Oceans modest increases in μmax result in higher δ15NNO3 due to enhanced local nutrient utilization, improving the agreement with reconstructions there. Models with moderately increased μmax fit both isotope data best, whereas large increases in nutrient utilization are inconsistent with nitrogen isotopes although they still fit the carbon isotopes reasonably well. The best fitting models reproduce major features of the glacial δ13CDIC, δ15N, and oxygen reconstructions while simulating increased Corg by 510–670 Pg compared with the preindustrial ocean. These results are consistent with the idea that the soft-tissue pump was more efficient during the LGM. Both circulation and biological nutrient utilization could contribute. However, these conclusions are preliminary given our idealized experiments, which do not consider changes in benthic denitrification and spatially inhomogenous changes in aeolian iron fluxes. The analysis illustrates interactions between the carbon and nitrogen cycles as well as the complementary constraints provided by their isotopes. Whereas carbon isotopes are sensitive to circulation changes and indicate well the three-dimensional Corg distribution, nitrogen isotopes are more sensitive to biological nutrient utilization

Inorganic carbon concentrating mechanisms in free-living and symbiotic dinoflagellates and chromerids

Photosynthetic dinoflagellates are ecologically and biogeochemically important in marine and freshwater environments. However, surprisingly little is known of how this group acquires inorganic carbon or how these diverse processes evolved. Consequently, how CO2 availability ultimately influences the success of dinoflagellates over space and time remains poorly resolved compared to other microalgal groups. Here we review the evidence. Photosynthetic core dinoflagellates have a Form II RuBisCO (replaced by Form IB or Form ID in derived dinoflagellates). The in vitro kinetics of the Form II RuBisCO from dinoflagellates are largely unknown, but dinoflagellates with Form II (and other) RuBisCOs have inorganic carbon concentrating mechanisms (CCMs), as indicated by in vivo internal inorganic C accumulation and affinity for external inorganic C. However, the location of the membrane(s) at which the essential active transport component(s) of the CCM occur(s) is (are) unresolved; isolation and characterization of functionally competent chloroplasts would help in this respect. Endosymbiotic Symbiodiniaceae (in Foraminifera, Acantharia, Radiolaria, Ciliata, Porifera, Acoela, Cnidaria, and Mollusca) obtain inorganic C by transport from seawater through host tissue. In corals this transport apparently provides an inorganic C concentration around the photobiont that obviates the need for photobiont CCM. This is not the case for tridacnid bivalves, medusae, or, possibly, Foraminifera. Overcoming these long-standing knowledge gaps relies on technical advances (e.g., the in vitro kinetics of Form II RuBisCO) that can functionally track the fate of inorganic C forms

Constraints on Earth system functioning at the Paleocene-Eocene Thermal Maximum from the marine silicon cycle

The Paleocene‐Eocene Thermal Maximum (PETM, ca. 56 Ma) is marked by a negative carbon isotope excursion (CIE) and increased global temperatures. The CIE is thought to result from the release of 13C‐depleted carbon, although the source(s) of carbon and triggers for its release, its rate of release, and the mechanisms by which the Earth system recovered are all debated. Many of the proposed mechanisms for the onset and recovery phases of the PETM make testable predictions about the marine silica cycle, making silicon isotope records a promising tool to address open questions about the PETM. We analyzed silicon isotope ratios (δ30Si) in radiolarian tests and sponge spicules from the Western North Atlantic (ODP Site 1051) across the PETM. Radiolarian δ30Si decreases by 0.6‰ from a background of 1‰ coeval with the CIE, while sponge δ30Si remains consistent at 0.2‰. Using a box model to test the Si cycle response to various scenarios, we find the data are best explained by a weak silicate weathering feedback, implying the recovery was mostly driven by nondiatom organic carbon burial, the other major long‐term carbon sink. We find no resolvable evidence for a volcanic trigger for carbon release, or for a change in regional oceanography. Better understanding of radiolarian Si isotope fractionation and more Si isotope records spanning the PETM are needed to confirm the global validity of these conclusions, but they highlight how the coupling between the silica and carbon cycles can be exploited to yield insight into the functioning of the Earth system